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ABSTRACT 

We present a detailed model of the discrete X-ray spectroscopic features expected from steady-state, 
low-density photoionized plasmas. We apply the Flexible Atomic Code (FAC) to calculate all of the nec- 
essary atomic data for the full range of ions relevant for the X-ray regime. These calculations have been 
incorporated into a simple model of a cone of ions irradiated by a point source located at its tip (now avail- 
able as the XSPEC model PHOTOION). For each ionic species in the cone, photoionization is balanced by 
recombination and ensuing radiative cascades, and photoexcitation of resonance transitions is balanced 
by radiative decay. This simple model is useful for diagnosing X-ray emission mechanisms, determining 
photoionization/photoexcitation/recombination rates, fitting temperatures and ionic emission measures, 
and probing geometrical properties (covering factor/column densities/radial filling factor/velocity dis- 
tributions) of absorbing/reemitting regions in photoionized plasmas. Such plasmas have akeady been 
observed in diverse astrophysical X-ray sources, including active galactic nuclei. X-ray binaries, cata- 
clysmic variables, and stellar winds of early-type stars, and may also provide a significant contribution to 
the X-ray spectra of gamma-ray-burst afterglows and the intergalactic medium. 

Subject headings: atomic data — atomic processes — fine: formation — plasmas — scattering — X-rays: 
general 



1. Introduction 

With the launch of the X-ray satellites Chandra (Weisskopf et al. 2002) and XMM-Newton (Jansen et al. 2001; 
den Herder et al. 2001), high-resolution X-ray spectroscopy of diverse astrophysical objects has now become routine. 
X-rays, due to their origin in extreme astrophysical environments as well as their significant penetrating ability, can 
potentially provide a wealth of information about these sources. Their penetration power, in particular, implies that 
X-ray emitting plasmas are often optically thin. In such cases, complicated radiative transfer is negligible, allowing 
for particularly simple astrophysical models. 

There are two main types of X-ray line-emitting plasmas: (1) Hot plasmas, which are mechanically heated 
through collisions and therefore have temperatures comparable to the observed line energies, and which produce 
X-ray line emission predominantly through radiative decay following electron impact excitation (hereafter, "colUsion- 
driven"), and (2) Photoionized plasmas, which are irradiated by a powerful external source and have subsequently 
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lower temperatures consistent with photon heating (including both photoelectric and Compton heating), and which 
produce X-ray line emission through recombination/radiative cascade following photoionization, and radiative decay 
following photoexcitation (hereafter, "radiation-driven"^ or "photoionized"). Hot plasmas are produced in the coro- 
nae of stars (including the Sun), in shock-heated environments (such as supemovae, cataclysmic variables, and stellar 
winds of early-type stars), and in the intracluster media of clusters of galaxies. Photoionized plasmas have recently 
been unambiguously confirmed in X-ray binaries (Liedahl & Paerels 1996, Cottam et al. 2001), AGN outflows (Sako 
et al. 2000, Kinkhabwala et al. 2002), cataclysmic variables (Mukai et al. 2003), and even in the stellar wind of the 
WR+O binary 7 Velorum (Dumm et al. 2002). Hybrid plasmas, in which both collision-driven and radiation-driven 
processes are important, can also exist. So far, hybrid X-ray-line-emitting plasmas have not yet been unambiguously 
detected in any astrophysical source, but they are expected to play an important role in accretion disks, for example, 
which are strongly irradiated (implying radiation-driven processes are important) but also have high densities (imply- 
ing collision-driven processes are important as well). Finally, we note that the mechanisms driving X-ray emission 
(either collision-driven or radiation-driven) in the two vastly different regimes of gamma-ray -burst (GRB) afterglows 
and the intergalactic medium (IGM) are still largely unknown. GRB afterglows may harbor hot plasmas produced 
through shocks and/or photoionized plasmas located in the ejecta or surrounding interstellar medium and powered 
by radiation from the burst and afterglow radiation. Similarly, the IGM may harbor hot plasmas from shocks due to 
structure formation and/or photoionized plasmas powered by radiation from active galactic nuclei (AGN) or starbursts. 

Calculations of X-ray line emission from hot, collision-driven plasmas have been discussed in great detail by 
several authors, including Raymond & Smith (1977), Mewe, Gronenschild, & van den Ooord (1985), Mewe, Lemen, 
& van den Oord (1986), Kaastra (1992), Liedahl, Osterheld, & Goldstein (1995), and Smith et al. (2001). Discussions 
of line formation in X-ray photoionized plasmas have also been presented by several authors, including Tarter, Tucker, 
& Salpeter (1969), Halpern & Grindlay (1980), Krolik, McKee, & Tarter (1981), Kallman & McCray (1982), Netzer 
(1993), Kallman et al. (1996), and Liedahl (1999); however, spectral predictions made in these works have not yet been 
compared in detail with real astrophysical X-ray spectra. In contrast, the assumptions and calculations we describe 
here for photoionized plasmas were made and tested in the process of analyzing high-resolution X-ray spectra from 
various astrophysical sources. 

Below, we present atomic structure and transition calculations using the pubUcly-available atomic code FAC 
(Gu 2002; Gu 2003) for all ionic transitions relevant for the X-ray regime. Atomic data values for some especially 
important transitions (in particular, more accurate wavelengths) and for some photoelectric edges are taken from other 
atomic databases, as will be described below. Using these data, we have fully characterized all significant line and edge 
absorption from the full range of included ions. We have also fully characterized the X-ray spectral emission features 
for radiative recombination forming the important H-Uke and He-Uke ions, as well as both radiative and dielectronic 
recombination forming L-shell ions. 

In order to attempt to explain real astrophysical spectra from photoionized plasmas, we have created a simple 
model of a cone of ions irradiated by a source located at its tip. Photoionization in the cone is balanced by recom- 
bination/radiative cascade and photoexcitation is balanced by radiative decay. This model has been incorporated into 
XSPEC (Amaud 1996) as the additive model PHOTOION^. For convenience, we have also created the abbreviated ad- 
ditive model PHSI for a single ion, and the multiplicative models MPABS (multi-ion absorption) and SIABS (single-ion 
absorption) for pure absorption studies. 



* "Radiation-driven" often refers to tlie force due to radiation pressure, but througliout tliis paper, tliis term is used only to refer to tlie ultimate 
power source for driving atomic transitions, subsuming both photoionization and photoexcitation. 

^Available at http://xmm.astro.columbia.edu/research.html 
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Our calculations invoke the following basic assumptions: (1) All ions are in their ground states. X-ray tran- 
sitions are generally much faster than the relevant rates for radiation-driven and collision-driven processes, making 
this typically a safe assumption. (2) ColUsional excitation and colhsional ionization are negligible. Due to the low 
temperatures and densities of photoionized plasmas, these processes are typically insignificant for driving X-ray tran- 
sitions. (3) We assume the entire medium is in a global steady state. (Note that we do not require the more restrictive 
assumption of local steady state conditions everywhere in the medium.) Furthermore, for the medium as a whole, 
we assume the total ionization rate equals the total recombination rate between all neighboring charge states. (4) All 
relative ion ratios are constant throughout the medium. Or, more approximately, all ion ratios are constant over regions 
with size comparable to the continuum optical-depth length scale (defined as the length over which the continuum is 
reduced by a few or more percent). For a mildly-absorbed medium, this length scale is equal to the entire medium, and 
this assumption is trivially valid. Interestingly, recent observations of the prototypical Seyfert 2 galaxy, NGC 1068, 
show that an inherent, relatively scale-free density distribution (over a few orders of magnitude) at each radius is more 
important than the radial distribution in setting the range in ionization parameter (Brinkman et al. 2002; Ogle et al. 
2003), suggesting that even in moderately or highly absorbed media this approximation may still be valid. (5) The 
irradiated medium can be modeled as a relatively narrow isotropic cone with the source of radiation located at its tip. 
This cone is completely specified by its opening angle, ionic column densities, representative ionic recombination 
temperatures, and a simple velocity structure, characterized by a global radial velocity shift and gaussian width, and 
a global transverse (perpendicular to the cone) velocity shift and gaussian width. (6) The emitting plasmas are opti- 
cally thin to their own emission. Because X-ray radiative cross-sections are generally weaker than longer wavelength 
cross-sections (e.g., in the UV or optical), a given medium remains optically thin to X-rays out to significantly higher 
column densities. Furthermore, in the context of our assumption of a narrow cone geometry, even with large radial 
ionic colunm densities and corresponding radial optical depth, the transverse optical depth is smaller than the radial 
optical depth by roughly the length-to-width ratio of the cone. Observationally, optical thinness has been demonstrated 
for the line emitting regions of NGC 1068 using a simple line conversion diagnostic (see Kinkhabwala et al. 2002 and 
§6.2 below). 

Our paper is organized as follows. Starting from the most general ionic rate equations in §2, we describe how 
the above assumptions are employed to simplify these equations for generic astrophysical plasmas in §3 and then 
for purely radiation-driven plasmas in §4. Here, we also describe our calculations of important diagnostics for H- 
and He-like ion emission, as well as the specifics of all of our atomic calculations with FAC. In order to explain 
real astrophysical spectra, we present our model of a cone of ions irradiated by a point source located at its tip in 
§5, obtaining the final expression for the rate equations, which serve as the basis for PHOTOION. In §6, we provide 
examples of the capabilities of PHOTOION. Finally, in §7, we discuss our results. 



2, Reaction Kinetics 

Presented below is a completely general analysis of the ionic rate equations for astrophysical plasmas (including 
all possible collision-driven and radiation-driven processes). This analysis is undertaken for two reasons. First, we 
are unaware of any similarly general analysis in the literature and hope to motivate studies along similar lines. And, 
second, this analysis most clearly brings out all the assumptions that are made in the formulation of our final rate 

equations for a radiation-driven medium. 

An extremely general form for the rate equation pertaining to population/depletion of ions nz,z,i of element Z 
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(atomic number) with z electrons in atomic level i and immersed in a generic gas or plasma is: 

k>\ j l<k<z j l<k<z j 

-n,,iY^ [R+ncC+{nC}Q+AY^j + ^n,j [R+n,C+ {nC}e+A]^;^. 

j j 

-nz.iY.Y.^^Ki'-^^.iY.Y^ii'^cjQy^i'^+j^ 

k>l j k>l j k>\ j 

Because only one element at a time is concerned, we have dropped the atomic number Z in the ion density subscripts 
on the RHS of the equation for conciseness, i.e., n^j = nz,z,j- The full time derivative on the LHS can be written as 
the usual convective derivative: dn/dt = dn/dt+y- Vn, which allows for the possibility of bulk motion (e.g., due to 
AGN outflow or stellar wind). Equation 1 includes all radiation-driven and collision-driven processes, i.e., all plau- 
sibly significant interactions with photons, electrons, ions, molecules, and dust grains and all possible excited-state 
decays. It therefore provides a good starting point for examination of most astrophysical plasmas, excluding only 
plasmas where coherent quantum effects are important, such as extremely-high-density plasmas (e.g., neutron star 
atmospheres) or coherent emission from masers or lasers (discussed further in the next paragraph). The notation is as 
follows: A (radiative decay rate), C ("collision-driven" coefficients), R ("radiation-driven" coefficients), D (autoion- 
ization rate), W (grain/molecule spontaneous disintegration rate), and parentheses (implying notation refers to all A, 
C, R, D, and W terms contained within), with the initial level as subscript and the final level as superscript for a 
particular transition between members of the ionic series (a la Einstein 1917). For example, C. / implies a transition 
from a z-hke ion (e.g., z = 2 for a He-like ion) in level i to a z'-hke ion in level j. {nC}Q refers to the set of prod- 
ucts of all collision-driven rate coefficients times their corresponding population densities n for interactions with ions, 
molecules, and dust grains (not including electrons). The right hand side, term-by-term, corresponds to (k+ l)-body 
recombination (z-kj) {z,i) (usually only two-body recombination is important, though three-body recombina- 
tion - i.e., two electrons and one ion - is significant at high densities), charge transfer (z-kJ) iz,i), ionization 
iz,i) — > iz-k,j), excitation/deexcitation/spontaneous decay iz,i) — > iz,j), excitation/deexcitation/spontaneous decay 
izj) (z,i), (k+ l)-body recombination (z,/) (z+kj), charge transfer (z,/) (z+k, j), ionization (z+k,j) — > iz,i), 
molecule/grain destruction q ^ q' + (z, 0^ molecule creation/grain adsorption q + (z, i) — > q' , and molecule/grain spon- 
taneous disintegration q ^ q' + (z, i). A similar rate equation to Eq. 1 but for molecules or grains is straightforward to 
write down as well. In light of the simpUstic approach and corresponding assumptions adopted in this paper, we do 
not investigate the coupling of Eq. 1 with any force laws, electrodynamic equations, or thermodynamic relations. 

The particular process of photoexcitation requires further elaboration. For photon-atom interactions near an 
atomic resonance (bound-bound), there are two possible regimes. One is the "classical" damped-oscillator regime 
(Wigner & Weisskopf 1930; Weisskopf 1931), in which it is possible to separately consider the quantum processes 
of absorption (photoexcitation) and reemission driven by spontaneous decay. This entire process can be thought of 
as "inelastic scattering." We have implicitly assumed this regime in our presentation of the ionic rate equations. The 
alternate regime is resonance scattering, in which these two processes cannot be separated out, and which always 
produces an outgoing photon with energy equal to the incident photon in the center-of-momentum frame ("elastic 
scattering"). To determine which is the valid regime requires consideration of the interaction timescale with the 
incident radiation field and the relative strength of the transition. Incident monochromatic light (formally infinite 
interaction time) yields energy-preserving resonance scattering for any transition strength. However, the continuous 
spectrum (short interaction time) or even absorbed power-law spectrum typical of photoionizing radiation and the 
predominance of absorption by relatively strong transitions appears to vaUdate the "classical" regime (e.g., Sakurai 



-5- 



1967). Observationally, the difference between the "classical" and resonance-scattering regimes is significant only 
in the details of the emission line profiles and in the expected degradation of higher-transition-energy photons to 
lower-transition-energy photons in the "classical" regime. 



2.1. Collision-driven Processes 

The C coefficients refer to all colUsion-driven processes, i.e., aU processes driven by colUsions of the particle of in- 
terest (usually an ion) with a secondary particle (usually a free electron). Collisional excitation/ionization/deexcitation 
and the processes of recombination, charge transfer, and molecule/grain creation/destruction can be compactly ex- 
pressed as: 

Ca=/ [a'']l(E)viE)fiE)dE, (2) 

where A and B, for example, could be (z,i) and (z'J). The integral is over the center-of-momentum energy, E, of 
the two particles. In the integrand, [cr'']^(£') is the collisional cross section for the transition A — > B of the particle 
under consideration (e.g., an ion) with the interacting particle (e.g., an electron), v(E) is the relative velocity of the two 
interacting particles, and, finally, /(£) is the normalized energy distribution. For isothermal electrons interacting with 
heavier, slower ions at the same temperature, E is the electron kinetic energy, v(£') ~ Ve(£') = ^2£'/me is the relative 
velocity, and f{E) is the electron energy distribution. Taking the MaxweU-Boltzmann distribution f{E;T) specified 
by temperature T for the electron distribution makes the final coefficient, C^{T), temperature dependent. 



2.2. Radiation-driven Processes 

The R coefficients refer to all radiation-driven processes, i.e., all processes driven by the local radiation field. 
These coefficients depend solely on the spectrum of the local radiation field and the given interaction cross section. 
Photoexcitation, photoionization, and stimulated emission and molecule/grain photodestruction can be expressed as: 

Rl{F{E,f)) = J [a'']l(E,r)F(E,r)dE, (3) 

where A and B again could simply be (z,0 and (z' The photon flux spectrum at position r is F(E,r) (e.g., with in 
units of photons cm"^ s^' keV~'; throughout this paper all luminosities and fluxes will be expressed in photon units, 
not the usual energy units). The photon flux spectrum may include radiation from a source external to the emission 
region as well as any other sources of radiation produced in the photoionized plasma itself. The total radiative cross 
section for transition A — » B is given by [i7^]^(£',f), which is already convolved with the local velocity distribution 
function at r of the particle population under consideration. Most of the calculational complexity of generic plasmas 
is reflected in F(E,r), which, through radiative transfer, is coupled to the density and velocity distributions (thermal 
or non-thermal) of all interacting particles. However, for the special case of optically thin plasmas, radiative transfer 
is neghgible, simpUfying calculations enormously (as is demonstrated below). 



3. Further Simplifications for Typical Astrophysical Plasmas 

We now make two assumptions valid for many typical astrophysical plasmas: (1) The only important collisional 
interactions involve an ion with a single free electron, and (2) Radiation-driven and coUision-driven ionizations remove 
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only one electron (this does not include subsequent electron ejection through autoionization, which is retained in the 
D coefficients). This gives for the ground-state equation: 



dt 



j>g 



j k>\ j 

and for the excited levels {i > g): 



^=j:^^rjn.c^rj- E E«.<'-«-E^:;/+E«^.i^S+E«-[^+«eC]:;; 

j l<k<z j j<i j>i j 



+E"^.^[^+"^'^]^j-"^.'E +E"^-i.^[^+"«<^]ztij+EE"^^*.^-^^47- 

i i j k>i j 

To predict the emergent spectrum, we must integrate Eqs. 4 and 5 over the entire emission region volume. For 
this purpose, we assume a global steady state, implying that all volume integrals / ^ c/V are identically zero (i.e., 
we balance ionic population/depletion and ionic inflow/outflow for the entire volume under consideration). A global 
steady state is assumed rather than a local steady state, because, for example, even in a steady outflow, the local rates of 
photoionization and recombination may not necessarily balance, with photoionizations occuring nearer to the source 
and recombinations occuring at larger distances. Even the assumption of a global steady state breaks down, of course, 
for highly-variable sources. However, under the global- steady- state assumption, integration of Eq. 4 over volume 
yields (with N = J ndV): 

i>g k>l j ^ j j 

-n,^i,,[/?+neC];f^_^ + «,.,^«eCf;'^-«,,,^ [R+n^CfAdV (6) 

j ' 

and integration of Eq. 5 over volume yields 

^.'(E4/+ E E^::^'^)-E^..'^:;^EE^^M^:tM= /(E«-'..«eC-^^ 

^ i<i \<k<z j ^ j>i k>\ i '' ^ j 

+Y,nzj [R+n,C] + E"zJ [R+ncC] ^'j + Y,n,^,j [R+n,C] . - n,,i [R+^A 'z!) ^V. (7) 
j J j j 

These equations are perfectly general for realistic collision- and/or radiation-driven plasmas. But they are still suf- 
ficiently complicated (due to coupling of these equations to each other and to other such equations) to bar direct 
integration. 



4. Radiation-Driven Emission 

In this section we show how further assumptions can be employed to decouple Eqs. 6 and 7 from each other and 
from all other such equations, leaving a simple analytic solution for derivation of the spectrum. But first we review the 
atomic processes relevant for radiation-driven plasmas. 
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4.1. Relevant Atomic Processes 

For radiation-driven plasmas, there are only a few processes which contribute significantly. These are summarized 
in Fig. 1 and explained below. 

The most important colUsion-driven process is recombination, including both radiative and dielectronic recom- 
bination. We presently discuss only radiative recombination. Using the principle of detailed balance embodied in the 
Milne relation (Milne 1924), the radiative recombination cross section for electrons at temperature T can be related to 
the photoionization cross section: 

-?l.(^e) = 7^^^^-:"'^(^th + £e), (8) 

where Eg is the electron kinetic energy, Eth is the threshold photonionizing photon energy, E = Eth+E^ is the photon 
energy, and the g are the ionic level degeneracies. Combining Eq. 8 with Eq. 2 allows for calculation of the radiative 
recombination coefficients. 

The two most important radiation-driven processes are photoexcitation and photoionization. For the photoexci- 
tation cross section, convolution of a line profile with the velocity distribution simply yields the Voigt profile (Voigt 
1912): 

7re2 1 1 

where Aud = v2 -)f and H(x,y) is the Voigt function with x = ^^^^^ (where A denotes the sum of all radiative 
and autoionizing decay rates from level i) and y = The Voigt profile is used for both absorption and reenussion 
line profiles (though with the possibility of different velocity distributions due to differing bulk motion distributions). 
The photoionization cross section a^~j^'^(E) is calculated exphcitly by FAC. We convolve the photoionization cross 
section with the same used above for the photoexcitation cross section. The total cross section for photoionization 
and photoexcitation is taken to be a sum of the individual cross sections. For photoexcitation, this simply means that 
we consider each transition (oscillator) decoupled from all other transitions (oscillators): A more accurate treatment 
would use the full Kramers-Heisenberg-like cross section with appropriate radiation damping (e.g., Sakurai 1967), 
but this level of accuracy is rarely, if ever, needed for astrophysical sources. The continuous and discrete line cross 
sections can therefore be taken one by one in Eq. 3 to calculate all of the necessary photoionization/photoexcitation 
coefficients. Finally, we note that photoexcitation and photoionization (as well as dielectronic recombination) can 
sometimes place electrons in levels capable of autoionization, as is illustrated in Fig. 1. 



4.2. Radiation-driven Rate Equations 

Returning to the ionic rate equations given in Eqs. 6 and 7, we now make two modifications. First, introducing 
a powerful simplification, we assume that all ionic excited levels are short-lived compared to the relevant radiation- 
driven and colhsion-driven excitation/ionization timescales. Second, in fine with purely radiation-driven plasmas, we 
assume that all colhsional excitations and ionizations are neghgible. We therefore obtain the following for Eq. 6: 

\ j ^ k>i j ^ j ^ 

and for Eq. 7: 



iX^^ E E<')-E^.^^^;^EE^-^.^^:4j= /("-i.^"'=c^i,^^ 

^ i<i \<k<z j ' j>i k>l j •' ^ ' 
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These equations have a few interesting properties, which we now briefly mention. In Eq. 10, the LHS contains only 
those transitions which connect charge state z to charge states with fewer electrons, and the RHS contains only those 
transitions which connect charge state z to charge states with more electrons. In Eq. II, all integrated expressions 
(involving only total ion numbers AO have been suggestively placed on the LHS with integral expressions placed 
on the RHS. The A,/?, and C coefficients in Eqs. 10 and 11 couple the same or nearest-neighbor charge states. The 
exception to this rule comes from the D terms (representing autoionization), which, through multiple electron ejection, 
can connect more distant charge states. 

Though most X-ray transitions occur very rapidly, allowing for our neglect of all excited-level populations, some 
rare long-lived levels do exist. The most important example is the especially long-lived 2^51 level in He-like ions (see 
§4.4 for more details). The inclusion of collision- or radiation-driven transitions out of these long-lived levels in this 
formahsm requires a much more sophisticated approach starting from Eq. 7. 



4.3. Radiation-Driven Emission from H-Like and He-Like Ions 

We now apply Eqs. 10 and II to the extremely important and relatively simple cases of H-like and He-like ions. 
Under the assumptions taken so far, autoionization is unimportant in He-like ions and onto H-like ions for the low 
temperatures characteristic of photoionized plasmas (and, of course, autoionization is trivially non-existent in H-like 
ions). We note that autoionizations, however, may be significant for driving ionizations onto He-like ions. With these 
considerations in mind, we can now simplify the rate equations for all charge states of interest (bare, H-like, and 
He-like). 

The absence of autoionization implies the only routes between bare, H-like, and He-like involve nearest neighbors 
(bare<-^H-like and H-like<-^He-like). Explicitly, for bare<-^H-like, taking z = in Eq. 10 yields: 

j ni,gRlg dV = j no^n^cl'j dV. (12) 

With a little more effort, a similar equation can also be derived for H-like^He-like transitions, which together with 
Eq. 12, can both be compactly expressed as follows (with z = I for bare<-^H-Uke and z = 2 for H-like-f-^He-like): 

J n,,gJ2Kg'' = / «z-l.««eCtotal dV (13) 

where Qotai = EjCz-i,g- 

Now turning to Eq. 1 1, for H-like ions (z = I electron), we obtain the following: 

^uY.^i-12^'-j'^J = J [non,C'/+ni,,R\:;^^dV. (14) 

Similarly, for He-like ions (z = 2 electrons) we obtain 

j<i j>i k>i j •' ^ ^ 

In both cases, we omit the irrelevant final term in Eq. II. The last term on the right-hand side of Eq. 15 gives 
the contribution to the particular He-like excited level i from autoionization of excited lower-ionization- state ions. 
Autoionizations ending on He-like excited levels are rare compared with those ending on the ground state (due to the 
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significantly higher energy change for n/ — > Is compared to nl n'l' with n' > 1). For this reason, we will henceforth 
neglect this contribution. Due to the symmetry of the H-like and He-like rate equations, we can now express Eq. 1 1 
as: 

with (with z = 1 for H-hke and z = 2 for He- like). Assuming that the temperature range is small enough that the 
branching ratios for recombination onto level i, C^l\ ^/Ctotai. are roughly temperature-independent, we can use Eq. 13 
to reexpress Eq. 16 as: 

}<i j>i J J 

Note that the RHS is purely dependent on the total photoexcitation/photoionization rates out of the ground state. For 
a given model of n^-^g and all R, the RHS of Eq. 17 predicts the individual (which are direct observables in the 
optically-thin limit). Also, since the source terms appear linearly on the right-hand side for photoionization (first 
term) with corresponding inverse process of radiative recombination/radiative cascade ("REC"), and photoexcitation 
(second term) with corresponding inverse process of radiative decay ("DEC") in Eq. 17, we can solve separately for 
each, obtaining Ni = Wf^*^ + Af^<^. 



4.4. Radiative Recombination in H-like and He-like Ions 

H-like and He-like ions dominate the reemission spectra of photoionized plasmas. Line emission in these species 
is also particularly simple. We therefore now discuss their important emission spectra in detail. 

The recombination spectrum expected for a total ionic recombination rate onto an ion with z - 1 electrons (ex- 
pressed as A^_j), is determined by the following equations: 

^.^E^::/-E^.^S=^^-i5^ (18) 

i j "^total 

with 

A^_l = y^nz-i,gneCtotal^^V~CtotalEMz_i, (19) 

where N^j denotes the number of ions z in level j, J^fj is the radiative decay rate for ion z and transition j ^ g 
{g denotes the ground state), C^l'j ^ is the recombination coefficient describing recombination of a free electron onto 
ground-state ion z- 1 creating ion z with electron configuration and Qotai = X^i^'fl'i g- For simplicity, we assume that 
recombinations occur at a single representative temperature (hence the ~ in Eq. 19). The emission measure EM^-i in 
Eq. 19 is simply defined as EM^-i = / nz-i,gnc dV. These equations are valid for both H-Uke and He-like ions (see 
Appendix). The set of equations represented by Eq. 18 provide a complete, soluble system of equations for all levels 
i > g. We note that the total population of ground-state ions is unobservable from the plasma fine emission spectrum 
alone; the best that can be done is to determine the emission measure EM^-i . 

Since we are interested in predicting line fluxes and ratios, we introduce the line luminosity (photons s~'): 



(20) 
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where Ij^i = NjA^^j/ A^_j is the dimensionless line luminosity coefficient. Similarly, the recombination luminosity 
onto level i is 

Lf^i = A:_,^ = Al,lf^,, (21) 

total 

where //^, = C^l\ ^/Qotai is the recombination luminosity coefficient (The subscript / stands for the initial state with 
a "free" electron). 

In Figs. 2 and 3, we plot the Une coefficients for the principal lines of H-like and He-like ions, respectively. Also 
shown in Figs. 4 and 5 are the corresponding Ctotal as a function of temperature. Together, these plots allow for the 
determination of the total recombination rates as well as ionic emission measures for a given photoionized plasma 
spectrum. In addition, we provide the standard He-like triplet ratios R = f/i and G = (;+ f)/r (Gabriel & Jordan 1969) 
for pure radiative recombination as a function of temperature in Fig. 6. We stress that since the He-like curves do 
not include the additional contributions from dielectronic recombination and collisional excitation present at the high 
temperature end, direct comparison with data for these high temperatures should be avoided. Explanation of exactly 
how these calculations were carried out is provided below in §4.6. 

We again stress the neglect of the well-observed process of forbidden-to-intercombination line conversion in the 
He-like triplets through excitation of the long-lived 2 ^Si level up to the 2 multiplet (Gabriel & Jordan 1969). 
Though we do not take this into account here (and it is not included in PHOTOION), it is simple enough to determine 
the effects on the He-hke triplet ratio R = f/i of electron colhsional excitation which is dependent on and (e.g., 
Porquet & Dubau 2000) and/or of UV photoexcitation which is dependent only on F(Eo) (e.g., Kahn et al. 2001). 



4.5. Radiation-Driven Emission from Other Ions 

Li-hke and lower-ionization- state ions are in general more complicated for a number of reasons. Excitations and 
ionizations can now lead to excited levels capable of autoionization of single and, importantly, multiple electrons. This 
implies that the strict equivalence between the single-electron removal rate (either through photoionization or single- 
electron autoioinization following photoexcitation) and the recombination rate between neighboring charge states 
are no longer necessarily equal. However, autoionizations removing more than one electron are typically negligible 
compared to photoionizations and single-electron autoionization (though this statement deserves further investigation). 
Therefore, this equivalence is likely still a very good approximation, and can be expressed as follows: 

= J «z-i.««eCtotai dV = J n,,,(Y,Kl' + T.^%f?,)jdV^ (22) 

where Qotai is the total recombination rate coefficient (now for both radiative and dielectronic recombination) and f^j 
denotes the fraction of decays from z, j that lead to electron removal via autoionization (hence the D superscript). 

We now focus on Eq. 11, which describes the line emission. As a further approximation, we drop the last term 
on the RHS of Eq. 11, which gives all autoionizations from lower charge states landing on the excited level z, i. This 
allows us to consider transitions only within the same ion. As for Eq. 17, we assume that the temperature range is 
small enough so that C^l'j / Qotai is independent of temperature. Then, from Eqs. 1 1 and 22, we obtain: 



^ j<i \<k<z j ' j>i tot J \ J J / 

+ / "z,sRlidV + 1 n,^,,,R%,JV. (23) 
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The linearity of Eq. 23 allows us to separate out the respective contributions to all A^^ , from photoionization/autoionization 
(eventually producing recombinations, both radiative and dielectronic) and the fraction of photoexcitations that pro- 
duce radiative decays: A'^,, = N^f''^ +A^™^. For photoionization/autoionization which produce eventual recombinations, 
we have: 

j<i j>i ^'"'^ ^ j j ^ 

We have neglected the J2i<k<z'^j^z^i' terms, which for the low temperatures typical of photoionized plasmas do 
not usually contribute significantly to the recombination spectra. These terms are retained, however, in Eq. 25, where 
they are essential for excitations and ionizations resulting in levels capable of autoionization. The system of equations 
represented by Eq. 24 can be thought of as the exact solution of the line spectrum for a given recombination rate onto 
each level i of A^_jC^l'j g/Qouh where A^_j is the total recombination rate. Similarly, for photoexcitation, we obtain 

^ j<i \<k<z j ' j>i 

Here, for the first time on the LHS, we must integrate over both n, g and n^+i ,^, the latter represent all photoionizations 
that leave the ion in an excited level (for photoionization of H- and He-like ions, this is always the ground state). As a 
further simplification to Eq. 25, we can discard the last term on the LHS, which describes feeding of z, i by all higher 
excitation levels: 

Nfr(X^zh E E^?') - ['^^,s^:sdV+ fn,.uR%.JV. (26) 
This directly gives the fluorescence yield for each excited level without having to perform a matrix inversion. 



4.6. Atomic Data Calculation and Implementation 

In the foregoing, we have dehneated how a series of approximations can reduce solution of the spectrum of each 
ion to a seemingly simple matrix inversion. We now discuss how these matrix inversions are currently implemented in 
PHOTOION. 

Throughout, we use atomic data generated by the atomic code FAC (Gu 2002; Gu 2003). For recombination 
calculations onto all relevant ions (for which recombination produces X-ray features), we explicitly calculate all levels 
up to principal quantum number n < 25 and all El, E2, Ml, and M2 transitions n'l' nl with n,n' < 25. Exphcit 
calculation of all levels for n = 45 aUows for interpolation for levels n = 26-AA. Using the H-like approximation for 
the recombination cross-sections, we sum all recombinations onto levels n > 45 and then spread this contribution 
over n = 45. Solution of the relevant matrices (such as in Eqs. 17 and 18) can now be obtained using only levels with 
n < 45. We also include all important two-photon decays for both H-hke and He-like (Drake 1986); such decays are not 
calculated by FAC. Similar calculations were carried out to detemiine the contribution from radiative and dielectronic 
recombination forming L-shell ions for all ions of interest in the X-ray regime. Due to the large number of excited 
levels for L-shell ions, inversions of the relevant matrix (represented by Eq. 24) required further approximations 
described in Gu (2002) and Gu (2003). 

Photoexcitation to all n < 25 levels {n< 100 for H- and He-like ions) is calculated directly from the atomic code 
results for the osciUator strengths and radiative decay rates Aji using the standard Voigt profile for line absorp- 
tion/emission discussed in §4.1. (Caveat: Since we do not presently calculate photoexcitation to higher n levels — 
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ideally out to n ^ oo — there is an artificial discontinuity at the photoelectric edges.) For radiative decay following 
valence-shell photoexcitation, solving the matrix for all excited levels is unnecessary since the dominant decay is the 
direct decay back to the ground state. Currently, to simpHfy and speed up the code, we therefore assume that every 
photoexcitation always produces a transition back to the ground state (even for photoexcitations to levels capable of 
autoionization, see below). This approximation is reliable for all important valence-shell transitions at the few percent 
level for reemission in an optically thin medium. For the same reason, we also only use the radiative decay rate to 
the ground for parameter x in the Voigt function (Eq. 9). Velocity widths in many X-ray absorbers have recently been 
observed to be significantly larger than the expected thermal width, which itself is much larger than the natural line 
width, making this a safe approximation. 

To increase the accuracy for especially important transitions (currently only in the H- and He-like ions), we veri- 
fied and, if necessary, modified important photoelectric cross sections, oscillator strengths, and especially wavelengths 
with previously-pubUshed calculations and experimental results from Vemer et al. (1996); Verner, Verner, & Ferland 
(1996); and NIST^ 

Also, not yet included is the calculation of all possible autoionization rates following photoexcitation or photoion- 
ization to determine total resulting line fluorescence. Currently, all photoexcitations are assumed to decay directly back 
to the ground state, leading to a sometimes gross overprediction of fluorescence for some excited levels with espe- 
cially high autoionization rates. Oppositely, fluorescence resulting from photoionization is currently absent in the 
code. Prospects for inclusion of these rates are discussed in §7. 



5. Irradiated Cone Model 

In order to integrate the expressions on the right-hand sides of Eqs. 17 and 23, we must first introduce our model. 
The simplest version consists of an isotropic cone irradiated by a nuclear source with photon luminosity spectrum 
L{E) =AE~^ (e.g., in units photons s~^ keV~'). We characterize the cone by a solid angle fl, radial colunm density 
G, and radial gaussian velocity distribution ct"'^ (which broadens all radiative cross-sections) for the ground state of 
each particular ion (Fig. 7). We first consider the more general non-isotropic case: L{E,9), G{6), and cr™''(^), where 
= i9,(l>). 

Employing the above assumptions, we obtain for the R coefficients pertaining to a particular ion n^^ (with z' = Z 
for photoexcitation and z' = z - 1 for photoionization): 

R^f(F(E,rM = J [a«]:/(£,^-)(^e-(^.^.«"))j£, (27) 

where T{E,r,9) = j^^^^, z}^'^^z,gi^' ^^)'^,ziE ,9)dr' and a§^iE,9) denotes the total radiative cross section for all 
radiation-driven transitions from the ground state of a particular element A with z electrons, with all cross sections 
already implicitly convolved with a possibly angular-dependent velocity distribution along the cone, al'^'^{9). L{E,9) 
is the intrinsic nuclear luminosity (again with possible angular dependence). The exponential takes into account ab- 
sorption along the cone due to the ground-state number densities n^^g. 

We now integrate the R coefficients times the appropriate ionic density over volume to get the total photoexcitation 
(or photoionization) rate in the plasma: 

J n,,g(r, 9)R('J(FiE, r, 9))dV = J dE J ^ J dr n,,g{r, 9) "'^;^{E,9) L{E, 9)e-^^''''-'^K (28) 



*http://physics.mst.gov/cgi-bin/AtData/main_asd 
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There are two different assumptions which can be made at this point to simplify this integral. The first assumption 
is that only the particular ion of interest contributes to the opacity. This is fairly robust if column densities of all other 
ions are low enough that only discrete line absorption is present in the spectrum, since the opacity due to lines alone is 
usually negligible. Under this assumption, the opacity r is only dependent on the particular ion of interest; therefore, 
dr = driE, r, 6) = n^^g{r, 9)a^(E, 6)dr, and we can rewrite Eq. 28 as: 

j n,,g{r,0)Rij(F{E,r,e))dV = J dE J ^ dr' L(E,e)e-^' . (29) 

Upon integration over r', this simplifies to: 



n^X'JdV = J —G,(9) J dE ^ ^(is, e") ^E, 6) , (30) 



-t(£,( 

t(eJ) 

where is the ground-state column density for ions (with z electrons) of the particular element under consideration. 



Alternatively, a second possible assumption could be that relative ionic abundances are independent of radius. 
Starting once again from Eq. 28, we now simply multiply and divide by the sum in large parentheses in the following 
expression: 



f ^ , f f dfl f n, Jr,tl)\a''\ (E,U) ^ ^\ ^ 

n,,,(r,e)RlME,r,e))dV= dE / dr /l k A E «z'..,,(r,^)4,z(£,^) U(£,^)^-^'^'^' 

(31) 

Recall that dT(E,r, d) = Yl{z -} ^2,z.g(x, 0)(J§ -(E ,d)dr. Therefore, the sum in large parentheses in Eq. 31 times dr is just 
dr. Therefore, by our central assumption that all relative ionic abundances are independent of r, the fractional term 
in Eq. 31 is also independent of r and can be replaced by a fraction involving only the total integrated ionic column 
densities Gz,z'. 



/ n,Jr,e)RlL\F{E,r,e))dV= / — / dE ^-^1^^^ L(E,e) / e"" dr'. 



(32) 



Gz(e)[a%;J(E,e) 
E{z,z}Gz,z(^")<,(£,^l 

Upon integration over r, we obtain: 

n,X:JdV = j j dE [a^]l '^\E,9)LiE,9)—^^j-, (33) 

which, interestingly, is identical to Eq. 30. 

For the simplest model of an isotropic cone of material with solid angle extent Q characterized by angle- 
independent radial colunon density G, radial gaussian velocity distribution a^, and nuclear photon luminosity spectrum 
L{E), the identical expressions derived above (Eqs. 30 and 33) reduce to: 

nz,,Rl:^dV = -G,j dE [a%'^^iE)L{E)—^. (34) 



/ 



This equation serves as the basis of PHOTOION. 



6. Spectra for Various Absorption and Reemission Geometries Calculated with PHOTOION 



The XSPEC model PHOTOION allows the calculation of spectra arising from different absorption and reemission 
geometries. Below, we provide examples of some simple geometries possible with PHOTOION. 
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6.1. Absorption 

An example of a spectram due to pure absorption is given in Fig. 8a. The imprinting of absorption features on 
the intrinsic continuum is determined by e'^^^^ with 

T(E) = T(E,0) = ! nz,z,gir',e)al^{E,e)dr' =Y, j nz,zA'-' ,^)4,z(E, 0)dr' = ^ Gz,zal,(E,0). (35) 

" {Z,Z} {Z,z}-^° {Z,z} 

Here, a§^^iE,6) denotes the total radiative cross section (in the line-of-sight direction 6) for all radiation-driven tran- 
sitions from the ground state of each ion, with all cross-sections implicitly convolved with the velocity distribution 
characterized by width cr'^''. This model can also be calculated using the related multiplicative XSPEC model MPABS 
(multi-phase absorber), which allows for pure absorption studies through the independent variation of each ionic col- 
umn density. These models should be useful in particular for understanding the complex, multi-phase absorption 
observed in a number of Seyfert 1 galaxy spectra (e.g., Kaastra et al. 2000, Sako et al. 2001, Behar, Sako, & Kahn 
2001). 



6.2. Reemission 

Observations transverse to the cone yield line spectra resulting from reemission in the cone. We concentrate on 
the dominant reemission from H- and He-like ions for simplicity (the contributions due to lower-charge-state ions is 
straightforward to include). The final spectrum of transitions to the ground state of ail relevant K-shell ions is then 
simply the sum of all ionic contributions: 

Reemission(£) = J^N.jA'^'jl j^i(E) + Al_,-^r(E), (36) 

^total 

where A^_j is the total recombination rate calculated from Eq. 13 and Ij^iiE) is the appropriate Voigt Une profile 

(normalized to 1) using velocity distribution cr[^'"'* and i'(E) is the ground-state recombination profile (normalized to 
1) convolved with aJ,™^ In general, cr^™** may not necessarily be equal to al'^'^. Explicitly, and without velocity 
convolution, the recombination profile is simply: 

r(E) cx f(E-E^;kn)^==apiiE), (37) 

where f(E-Exb;kTe) is the Maxwellian distribution, evaluated at £-£th and specified by temperature kT^, and crpi(£) 
is the photoionization cross section, with cr^iiE) = for £ < Eth- The contribution from radiative and dielectronic 
recombination forming lower-ionization-state ions is also straightforward to include. Examples of reemission spectra 
for differing amounts of obscuration towards the intrinsic continuum are shown in Fig. 9). 

If the emission regions are moderately optically thick to their emission, then the reemitted spectra will be modi- 
fied. Higher-order transitions (n ^ 1 with n>2) will be degraded through multiple scattering to lower-order transitions 
(at sufficiently high optical depths, all higher order lines are converted to n = 2 ^ 1 transitions). This reprocessing 
can provide an estimate of the overall optical thickness of the medium. In Fig. 13 of Kinkhabwala et al. (2002), the 
diagnostic value of this effect in H-like and He-like line series is illustrated. The lines shown in that plot give the 
column density and velocity width in (assuming a roughly spherical medium) corresponding to 10% conversion of 
the Ly/3 photons to Lya and, similarly, 10% conversion of He/3 to Hea transitions for C through Mg (degradation 
percentage can be scaled arbitrarily in these plots). 
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6.3. Absorption Plus Reemission 

Even if the ionization cone is viewed in absorption, there is some amount of the reprocessed emission which will 
be reemitted in the observer's direction (see Fig. 7). The total spectrum can be written as SpectrumCE) = L(£')e~^^^^ + 
ReemissionCii), where ReemissionCf) = Decay(£') + Recombination(£'), with DecayC/i) giving the total contribution 
from radiative decay following photoexcitation and Recombination(£) giving the total contribution of both radiative 
and dielectronic recombination. Both DecayCE) and Reemission(£') impUcitiy contain the details of absorption in the 
cone, as shown below. 

In order to account for absorption of the reemitted spectrum in the cone we must make some further assumptions. 
First, we assume that the cone is sufficiently narrow so that the absorption factor pertaining to any region of the cone 
is just dependent on radius. We define two such absorption factors: To(£, r) = J^^ J^iz z} ^z.z.g(r) [ct*] ^ - g(E), repre- 
senting absorption from the irmer radius Vg of the cone to r, and Tf{E, r) = J^' z} ^z,z,g(T) [f*] ^ ^ J^^)^ representing 
absorption from r to the outer radius of the cone Vf. Clearly, the total opacity is just t{E) = ToiE,r) + TfiE, r) for any 
value of r. 

First, we derive a general expression for Decay(£'). For simplicity, we assume that photoexcitation corresponds 
to coherent resonant scattering (in contrast to statements made in §2). The relevant spectrum expressed as a simple 
integral over the cone can then be written as: 

Decay(£) = J] /'F(£,r)nz,.,g(r)^ [a*]^';'^(£)e-"/W^2^^jQ = 

^^LiE).-'''Y.T.y&^^ ^ 

{z,z} i 

ilL(£)e-<-'EEG-.[-1z;:^(^)- (38) 

{Z,z} j 



Accounting for absorption of the reemission contribution from recombination is more difficult due to the need 
to integrate over each photoelectric edge to determine the photoionization (and, therefore, recombination) rate. Here, 
we define Rz.ziE) to be the spectrum (properly normalized) of recombination (radiative and dielectronic) forming ion 
Z,z. Total recombinations are balanced with photoionizations (resulting either from direct photoionization or indirect 
ionization via autoionization following photoexcitation; though only the contribution from the former is shown below 
for simplicity). This yields: 

Recombination(£) = ^Rz,z{E) j i j FiE\r)nz^,gir)^[a'^]^/~^'\E')dE'je-^f'^'^r^drdn = 

{Z,z.} •' j '^"^ ^ 

^ ^ RzAE) I ( / L(£')e--(^'^'' ^ [a^ '^-l'\E')dE) ^-^f^^W,,,ir)dr. (39) 

Here, the integral over energy in parentheses is still a function of r through To(r,E), requiring a costly double integral. 
In order to speed up this calculation, we would Uke to spht this double integral into two sequential integrations over 
E and then r. This requires approximating the opacity term within the energy integral as some average opacity 
independent of radius, e"^"^'''^^ « Opacity(£). Employing once again the assumption that all ionic ratios are constant 
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throughout the cone, we obtain: 

RecombinationCE) =^Y1 ^^-^^^^ ( / LiE')OpacnyiE')Y^ [a'^f^l'^'\E')dE'j x (40) 
{z,z} j '''^ ^ 

I ''^''""'^ ^ n,,,,ir)4>,iE)dr (41) 

J Z^{z',j}«z',z,g('-)CTz,_,(ii)^^^ 

Recalling that dTf{r, E)= — X^^^' z} '^z' ,z,gi'')(^' z(E)dr, we can integrate this expression (e.g., see Eqs. 3 1-33), obtain- 
ing: 

Recombination(£) = ^ ^ Gz,zRz.z(E) "^^^ / L(£')Opacity(£') ^ [a^] ^^'~^'\E')dE' . (42) 
{z,j} "Tl. ) 7 ^. 

Useful lower and upper limits to the recombination contribution are obtained by taking Opacity(£') = e~^^^^ and 
Opacity(£') = 1, respectively. Spectra resulting from assumption of either limit can be calculated with PHOTOION 
(see Fig. 8). 



7. Discussion 

The models presented in this paper have already received experimental verification through X-ray spectra of 
multiple accretion-powered astrophysical sources. In particular, application of our model to the remarkable XMM- 
Newton spectrum of the brightest, prototypical Seyfert 2 galaxy, NGC 1068 (Kinkhabwala et al. 2002), provided the 
first quantitative test of a photoionization code (e.g., the gratifying confirmation of the ratio of intercombination line 
to RRC in O VII). In addition, the self-consistent inclusion of photoexcitation and resulting radiative decay allowed 
for discrimination between emission mechanisms in the spectrum of NGC 1068, providing a robust upper limit to any 
additional contribution from hot plasmas. Observations of NGC 1068 also provided the motivation for the seemingly 
ad hoc assumption taken in §5 that all ionic ratios are constant throughout the cone (Brinkman et al. 2002; Ogle 
et al. 2003), as was explained in §1. We expect that PHOTOION (and associated codes PHSI, MPABS, and SIABS) 
will continue to be of use in the description of astrophysical photoionized plasmas existing in a variety of sources, 
including AGN, X-ray binaries, cataclysmic variables, stellar winds of early-type stars, GRB afterglows, and the IGM. 

Lastly, as mentioned in §4.6, radiative and autoionizing decay rates (as well as the resulting fluorescence intensi- 
ties) for levels capable of autoionization have not yet been properly included. EspeciaUy important are the Fe K-shell 
transitions in Fe M-shell and L-shell ions. The number of atomic levels involved in these calculations is large, but not 
prohibitive. These represent the final atomic calculations necessary to "complete" PHOTOION, at least in terms of the 
assumptions taken in this paper. These rates should eventually be included in a future version of the code. 
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Fig. 1 . — Highly- simplified grotrian diagrams depicting all of the relevant atomic processes for photoionized plasmas. 
The top two diagrams depict photoexcitation and its inverse process of radiative decay (only the dominant direct 
decay path is shown, but other decay paths are possible). The next two diagrams depict photoionization and the 
inverse process of radiative recombination and radiative cascade. Another important inverse process to photoionization 
(especially for multi-electron ions) is the two-step process of dielectronic recombination, depicted in the next two 
diagrams. The last three diagrams depict photoexcitation up to a level capable of autoionization followed by either 
radiative decay or autoionization; these diagrams can be trivially modified to describe photoionization to levels capable 
of autoionization as well. 
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Fig. 2. — Line coefficients for the brightest features resulting from radiative recombination forming H-like ions. Ly a 
(black), Ly f3 (cyan), Ly 7 (green), Ly S (olive), and Ly e (orange) are shown, along with the ground-state RRC (red). 
Division of the observed line luminosity by the line coefficient gives the radiative recombination rate A^_[ (at an 
assumed temperature). Division of A:_i by Ctotai in Fig. 4 gives the ionic emission measure EM,_i = J 11^,-1. grigdV (see 
Eq. 19). 
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Fig. 3. — Line coefficients for the brightest features resulting from pure radiative recombination forming He-like ions 
in a low density plasma. Forbidden (black), intercombination (purple), resonance (blue), He/3 (cyan), He7 (green), HeJ 
(olive), and Hee (orange) are shown as well as the ground-state RRC (red). Division of the observed line luminosity 
by the line coefficient gives the radiative recombination rate A;_[. Division of A^_j by Ctotai in Fig. 5 gives the ionic 
emission measure EMj_i = j n^-i^gHf-dV (seeEq. 19). Caveat: At the low temperature end, only radiative recombination 
is important, but at the high temperature end {kT ^ Eq where is the resonance line energy), additional contributions 
to the line emission (not included) from dielectronic recombination and collisional excitation become important. Also, 
high densities or strong ambient UV fields can convert / i, as explained in §4.4. 
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Fig. 4. — The temperature dependence of Ctotai for radiative recombination forming H-like ions (see also Fig. 2). 
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Fig. 5. — The temperature dependence of Qotai for radiative recombination forming He-like ions (see also Fig. 3). 
Dielectronic recombination, which becomes important at the high temperature end {kT ~ Eo where Eq is the resonance 
line energy), has not been included. 
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Fig. 6. — He-like triplet ratios R = f /i and G = (i+ f)/r for pure radiative recombination for elements C, N, O, Ne, 
Mg, Al, Si, S, Ar, Ca, Fe, and Ni (from top to bottom). At the low temperature end, only radiative recombination is 
important, but at higher temperatures {kT ^ Eq where Eq is the resonance line energy), additional contributions (not 
included) from dielectronic recombination and coUisional excitation become important. 
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Absorption 




Fig. 7. — Simple cartoon of irradiated cone model (not to scale). The central nuclear component which irradiates 
the cone is shown as the black spot. The labels "Absorption" or "Reemission" denote the two main orientations 
of interest. The following model parameters are indicated: f = fl /A-tt (covering factor), A^'^^ (each individual ionic 
column density), a'f'^ (radial gaussian velocity width), ct[™^ (transverse gaussian velocity width), and Lx (irradiating 
spectrum). In the "Absorption" direction, in addition to features imprinted by line and edge absorption, reemission in 
the cone may also provide an additional contribution (see Fig. 8). Also, in the "Reemission" direction, the intrinsic 
continuum may or may not be obscured by an additional medium (see Fig. 9). Finally, we note that some sources may 
also exhibit a second cone (shown at the bottom of the figure). 
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Fig. 8. — Three examples of multi-phase "Absorption" spectra (see Fig. 7) calculated by PHOTOION: a) Pure absorp- 
tion, b) Absorption plus reemission (upper limit assuming Opacity(£') = 1 , see §6.3), and c) Absorption plus reemission 
(lower limit assuming Opacity(£') = e"'"'^^*, see §6.3). For all spectra, we take cr™'' = 400 km s~' and a power-law slope 
of r = 1.8. And, for spectra in panels b) and c), we further assume a covering fraction of / = il/Air = 0.15 for the 
reemission component. 
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Fig. 9. — "Reemission" analogue (see Fig. 7) to the same absorption model used in the top panel of Fig. 8. For all 
spectra, we take aj^^'' = aj,™*' = 400 km s~\ r= 1.8, and /= r2/47r = 0.15. Here the three panels correspond to differing 
amounts of obscuration along the observer's Une-of-sight to the intrinsic continuum: a) Pure reemission spectrum 
(intrinsic continuum is completely absorbed), b) Reemission spectrum plus moderately-absorbed intrinsic continuum 
(with neutral opacity towards intrinsic continuum equivalent to A^h = 5 .0 e22 cm~^), and c) Reemission plus unabsorbed 
intrinsic continuum. 
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